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INTRODUCTION 



I . 

This study was part of an ongoing project to produce an 
empirical heat budget parameterization for utilization in the 
Fleet Numerical Weather Central (FNWC) primitive equation 
prediction model. The project will eventually examine the 
heat balance of the earth-ocean system using representative 
mid-seasonal data from each of the four seasons of the year 
1973-74. This particular study covers 15-16 October 1973, 
using gridded FNWC temperature and humidity fields, at con- 
stant pressure levels, for this period. 

The parameters which influence the radiative calculations 
most strongly are the two cloud layers of fractional amounts 
CL(1) and CL(2), at 600 and 900 mb, respectively, at each 
grid point. The fractional amounts are specified using the 
large-scale formulations developed by Smagorinsky (1960) as 
also was done by Jenks (1974). The possibility of two cloud 
layers introduces substantial complexity to the radiative 
physics compared to the one-layer cloud model previously de- 
scribed by Martin (1972). The adaptation of the radiational 
model to FNWC gridpoint soundings is based on an unpublished 
manuscript (Martin, 1974). The Martin manuscript in turn 
has similarities to the multi-cloud heating packages of 
Arakawa (1972) and Rodgers (1967). 

Heat-balance computations were made at the ocean surface, 
at the tropopause, and for the atmospheric layer surface to 
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200 mb. A major purpose was to determine whether the cloud- 
cover parameterization previously adopted by Jenks (1974) 
lends itself to realistic budgetary computations (as compared 
to the results of others) for Oct. 15-16 over the Northern 
Hemispheric oceans. In studying the net heating effects at 
the ocean surface, the evaporative and sensible-heat trans- 
fers were modeled after Kesel and Winninghoff (1972), and 
Kaitala (1974), but with some modifications to be described 
later (Section V). 

To assess the effects of the cloud cover parameterization 
on the heat transfers of this study, gridpoint data along 
four meridians were selected from the FNWC analyses of 15-16 
October 1973. Three of the meridians are from the Pacific 
Ocean area and the fourth is from the Atlantic Ocean (see 
Fig. 2). Most of the gridpoint soundings sampled were in 
the Northern Hemisphere. 

The five sigma (a = p/tt) analysis levels of the FNWC 
( primitive equation model were adapted for data input in the 
vertical scale of this study. The k-levels (k=l, . . . , 10) , 
where k=10a, are used for the purpose of simplified refer- 
encing of the data levels in the equations derived and are 
shown in Fig. 1. The cloud layers, when they exist, are con- 
sidered to be in two layers, a high-level cloud layer from 
k=4 to k=6 and a low-level cloud deck from k=8 to k=9, with 
amounts CL(1) and CL(2), respectively. In pressure-scaling 
the water vapor and C0 2 absorber masses a realistic method, 
after. Smith (1966), was used to extrapolate water vapor 
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-0.7 



0= 0.8 



( 8 , 10 ) 



T.q- 



0.9 



-i- 10 0=1.0 

Figure 1. Five-layer radiative sounding used in this study. 
Levels are identified by their values on the k-seale, while 
layers are identified by their level boundary indices in 
parentheses, o.g. (8,10). Pressure-scaled water vapor and 
C0 2 mass increments U and C, respectively, are integrated 
with respect to the surface and introduced at even levels 
while the mixing ratio, q, is formulated at odd levels. The 
temperature, T, is expressed for all levels. 
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mixing ratio values above the level of k=4 , where radiosonde 
moisture data is normally not available. The net flux of 
terrestrial radiation computations make use of empirical ab- 
sorptive, ties of water vapor and C0 2 (Sasamori, 1968), which 
are incorporated in the Martin (1974) manuscript. The solar 
radiation package in Martin (1974) has been modeled after 
Arakawa (1972), except that only two cloud layers are con- 
sidered possible, as contrasted with three in the Arakawa 
version . 
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II. DATA PREPARATION 



The data used in this study was in the form of gridpoint 
soundings chosen along four meridians of FNWC 63-by-63 grid- 
mesh temperature and humidity data fields to insure minimum 
over land data. The three meridians (and their respective 
number of soundings) selected over the Pacific Ocean were 
located at 125°W - (25 soundings), 170°W - (25 soundings), 
and 145°E - (17 soundings). The Atlantic Ocean meridian was 
35°W - (26 soundings). Fig. 2 shows how the meridians were 
located on the FNWC polar stereographic grid. This procedure 
of choosing data lines along diagonals or along J=32 (through 
the pole-point) of the FNWC data-fields obviates the neces- 
sity of interpolating spatially to gridpoints. Data along 
line 3 in the Pacific were not extended southward of (9,55) 
because they fell over land masses (New Guinea and Northern 
Australia) where the surface temperatures were too high to 
be representative of the oceanic values. 

The vertical soundings on these data lines were taken 
from the original FNWC 63-by-63 gridpoint analysis of tem- 
perature at the surface and at nine other standard levels up 
to 100 millibars for 0000GMT , 16 October 1973 for the Pacific 
Ocean data cases and 12 hours earlier for the Atlantic Ocean 
data cases. In addition to these standard level values, a 
T , value was listed which is the FNWC calculated temperature 

X 

(see Kesel and Winninghoff, 1972) for the top of the turbu- 
lent layer of constant sensible heat flux. This level, 
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( 1 . 63 ) 




Figure 2. FNWC polar stereographic grid and meridians (lines 
1, 2, 3, and 4) selected for study. The longitudes A are 
shown for each meridian as well as the extent considered of 
each meridian. 
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which is approximately 20 meters above the surface, was des- 
ignated by a code label of 9999 under the column labeled 
pressure (Table 1(a)). The six moisture values of this table 
were given by the vapor pressure at sea level, and by dew- 
point depression for levels from 925 mb to 400 mb. In addi- 
tion, the parameter e was provided by FNWC as the vapor 

X 

pressure at the top of the constant flux layer at the level 
indicated by the code 9999 (Table 1(a)). Table 1(a) gives 
an example of the original sounding for gridpoint (1,1) of 
the FNWC map for the given period. 

In order to perform radiative calculations, it is neces- 
sary to have water vapor and C0 2 absorber masses and cloud 
amounts within the required k-level boundaries (Fig. 1). To 
do this, a systematic method was applied to arrive at the 
temperature and mixing ratio distributions (Table 1(b)) as 
detailed in the following subsections. It should be noted 
that all soundings start at sea level, but as a reasonable 
approximation surface pressure was considered to be tt^IOOO mb. 
Therefore the k-levels of Fig. 1 become 1000., 900., 800., 
...., 200., 100., 0.0 mb, respectively. 
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Tabic 1(a) . Example of an original FNWC gridpoint sounding 
for point (1,1). The humidity parameters for the surface 
and the tope of the constant flux layer (9999) are vapor pres- 
sures and those between 925,..., 400 mb are dew-point depres- 
sion . 



Pressure (mb) 


T(°C) 


Humidity Parameters 


surface 


25.5 


25.6 mb 


925 


19.3 


3 . 6° C 


850 


15.5 


3 . 2° C 


700 


6.5 


5 . 4° C 


500 


-10.2 


o 

o 

lO 


400 


-21.9 


2. 1°C 


300 


00 

CD 

00 

1 




250 


1 

OI 

00 




200 


-56.0 




150 


-67.3 




100 


CD 

O 

00 

1 




9999* 


T = 24.7 

X 


e = 26.3 mb 

X 



Table 1(b). 


Example 


of the corresponding 


radiative 


sounding 


with temperature and 


mixing ratio 


listed ; 


at k-levels 


• 


Pressure 


(mb) 




T(°K) 


Mixing Ratio 


(g/kg) 


1000 






298.7 




15.92 




900 






290.9 




11.45 




800 






285.9 








700 






279.7 




5.95 




600 






272.4 








500 






263.0 




2.31 




400 






251.3 




1.42 




300 






236.4 




0.56 




200 






217.2 








100 






192.6 




0.03 




0 






192.6 




0.00 




* Values at 


top of 


constant flux 


layer as 


used in FNWC heat- 


ing package 


(Kesel 


and Winninghof f , 1972) 


. 
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A . TEMPERATURE TERMS 



The temperatures were listed at each mandatory level be- 
tween the surface and 100 mb. The temperature T 10 was set 
equal to the listed surface air temperature. At the top of 
the atmosphere, the temperature was assumed to be isothermal 
from 100 mb to P=0 mb, with the mixing ratio approaching 
zero at P=0 mb. The radiative sounding temperatures for the 
remaining k-levels were obtained from either their corre- 
sponding listed temperature or by a three-point Lagrangian 
interpolation scheme of the form: 

T k = T o(( P k - Pl)(P k“ P2))/((Po_Pl)(Po_Pl)) 

+ T 1 ((P k -P 0 )(P k -P 2 ))/((P 1 -P 0 )(Pi-P 2 )) (2-1) 

+ T 2 ((P k -P 0 )(P k -Pi ))/((P 2 -Po XP 2 -P 1 )) 

Here k is the pressure level desired for T k interpolation, 

P 0 lies above P k , and Pi, P 2 lie below level k. The param- 
eters T 0 , T) , and T 2 are temperatures at the same levels as 
P 0 , Pj and P 2 , respectively. As an example, if k=9, Pi and 
P 2 were taken as 925 mb and 1000 mb, respectively. 

The T^ temperature was not utilized in the interpolation 
scheme described above, but was carried along as a twelfth 
temperature for later use in the heat budget computations. 

B. MOISTURE TREATMENT IN THE SOUNDING 

The moisture parameters from the original FNWC soundings 
(Table 1(a)) were converted into mixing ratios at each of 
the original sounding levels. At k=10 the surface vapor 
pressure, e a ^ r » was used to calculate the mixing ratio as 
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follows : 



Qi o o o : 


* .62197 e . /1000 (2-2) 

9.ir 


where e . 

air 


is in millibars. To obtain a mixing ratio for 


the top of 


the constant flux layer, Equation (2-3) was used 


to compute 


q from e 

^x X 



q - (e /1000) .62197 (2-3) 

X X 

where e is the vapor pressure at the top of the constant 

X 

flux layer in mb. 

The mixing ratio at level P is calculated at the remain- 
ing original data levels by the use of Equation (2-4) which 
is taken from Fleagle and Businger (p. 63, 1963); 



q(P) = 


621 97 e (A " B/T) -B(T-T ) 

- p exp[ TT ^ D ] (2-4) 


A 


= 21.656 


B 


= 5418. 0°K 


T 


= temperature °K at level P 



Tp = dew point (°K) at level P 

The calculated q-values were then interpolated to the k-lev- 
els using an equation similar to Equation (2-1). An example 
of the interpolation results for gridpoint (1,1) is shown in 
Table 1(b). 

Since accurate humidity information is not generally 
available for the upper troposphere, a logarithmic extrapola 
tion formula was used to find q-values at k=l, 2, and 3. 



sin = 

qs 


<500>* C2 - 5) 
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This formula as described by Jenks (1974) is similar to that 



developed by Smith (1966). It is solved, using a "best-fit" 



eter X. The six q- values of the original soundings corre- 
sponding to P = 1000, 925, 850, 700, 500, and 400 mb are 
used to determine X, which is then used in Equation (2-5) to 
extrapolate the q values to P = 300, 200, and 100 mb. 

A correlation coefficient was statistically computed in 
connection with each X-profile of the form: 



The correlation coefficient enabled the monitoring of the 
usefulness of the X-profile technique and was generally 
found to be in the range of 0.95 to 0.99. 

C. PRESSURE-SCALED ABSORBER MASSES 

Having computed the mixing ratios at each odd k-level, 
it is now possible to compute the pressure-scaled water va- 
por absorber masses in each layer. In Figure 1, these water 
vapor layer-masses were denoted by U, however here we will 
use notation which is more representative of the particular 
layer under consideration. For example, the pressure-scaled 
water vapor mass in the lowest layer (8,10) is given by 



a P 

between variables y = log and x = log for the param- 

Q 5 bUU 





( 2 - 6 ) 
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where the Moller and Raschke (1964) pressure-scaling exponent 
0.72 has been employed. For the layer (6,8), an analog to 
Equation (2-6) was used by replacing the parameters with the 
subscript 9 with their counterpart parameters with the sub- 
script 7. A similar procedure was followed for the remain- 
ing layers (4,6), (2,4), (0,2). In each of the five layers 
AP = 200 mb. 

The algorithm for computing the integrated water vapor 
mass above the earth's surface is as follows: 



U(10) 


= 0.0 




U( 8) 


= U( 10) 


+ U( 8 , 10 ) 


U( 6 ) 


= u(8) 


+ U( 6 , 8 ) 


U(4) 


= U(6) 


+ U( 4 , 6 ) 


U( 2 ) 


= U(4 ) 


+ U( 2 , 4 ) 


U(l) 


= U(2) 


+ 0 . 5U ( 0 , 2 ) 


U(0) 


= U(2) 


+ U( 0 , 2 ) 



A similar computation scheme is used for the carbon di- 
oxide scaled mass. Here, however it is customary to repre- 
sent the pressure-scaled mass in terms of the normal 
temperature and pressure (N.T.P.) volume over a cm 2 of the 
earth's surface. Thus for the lower layer: 

C(8, 10) = 3.14 x 10-'(||-)( TM §* 1I )”'’ 2 (2-8) 

s 

The fractional amount 3.14 x 10 -1 * is the constant mixing ra- 

AP 

tio of CO 2 by volume, and is the N.T.P. thickness of the 

• gP s 

layer AP = Pio-Pe- The quantity p is constant density in a 

s 

standard homogeneous atmosphere of depth II where T =273.16 K 
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everywhere and P = 1013.25 mb. The parameter H, also gen- 
erally known as the scale height of the homogeneous atmos- 
phere, may be further defined as: 

RT 

H = = 7.995 x 10 s cm 

g 

Therefore Equation (2-8) may be rewritten as 

C(8,10) - 3.14 x 10- (2-9) 

All other layer thicknesses of C0 2 are obtained in a similar 
manner to Equation (2-9) by substituting the appropriate 
subscripted indices 7,5,...,1 on the pressure P. 

The algorithm for computing the integrated carbon diox- 
ide mass above the earth's surface is as follows: 



C(10) 


= 0.0 




C(8) 


= C( 10) 


+ C( 8 , 10) 


C(6) 


= C( 8 ) 


+ C( 6 , 8 ) 


C(4 ) 


= C( 6 ) 


+ C(4,6) 


C( 2 ) 


= C(4 ) 


+ C(2,4) 


C(l) 


= C( 2 ) 


+ 0. 5C(0, 2) 


C(0) 


= C(2) 


+ C( 0 , 2 ) 



D. CLOUD PARAMETERIZATION 

In order to determine the fractional cloud cover amounts 
in the layers (4,6) and (8,9), the saturated vapor pressure 
must be calculated at levels k=5 and k=9, respectively by 
the following equations: 



q s (500) 



.62197 e (A B / Ts) 
500 



(2-lla) 
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q s (900) - 



.62197 e (A B / T9 * 



900 



(2-llb) 



where A = 21.656 and B = 5418. 0°K as before in Eq. (2-4). 

Knowing the saturated vapor pressure at levels k=5 and 
k=9 enables the determination of the relative humidity at 
those levels: 



R H ( 5 ^ = Q ( 500 ) 
R ' H ' (5) q s (500) 

R H f Q ^ = Q C 900 ) 
R ' H,(9) q s (900) 



(2-12a) 

(2-12b) 



The relative humidities were then substituted into 
Smagorinsky ’ s (1960) formulation for fractional cloud amounts 



CL( 1 ) = 2.0 (RH( 5 ) ) - 0.7 



CL(2 ) = 3.33 ( RH ( 9 ) ) - 2.0 



(2-13a) 

(2-13b) 



These formulations were employed by Jenks (1974) and are 
also to be tested in this study in which a different season- 
al sample is considered. 

These fractional cloud amounts are a function of large 
scale influences only and do not consider small scale con- 
vective activity, seasonal conditions or latitude. The use 
of such clouds is for radiational study only. 
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III. TERRESTRIAL RADIATION 



A. THEORETICAL AND EMPIRICAL BASIS 

Empirical formulas were developed by Sasamori (1968) for 
computing emissivities for long-wave flux calculations asso- 
ciated with the NCAR General Circulation Model. These em- 
pirical emissivity formulas can be related to the values 
built into Yamamoto's (1952) well known Radiation Chart, 
which was developed from radiative transfer theory. This 
chart has proved to be very accurate for computational pur- 
poses, and was used by Jenks (1974) as a schematic guide for 
integration of the radiative transfers through the various 
layers of interest in the soundings discussed in Section II. 

The necessary emissivity formulas of Sasamori for utili- 
zation in the FNWC heating package are those which will yield 
expressions for: 

F 10 * = Net IR flux at earth surface 
F 6 * = Net IR flux at level k=6 

F 2 * = Net IR flux at level k=2 

F610 = Net IR flux divergence in the layer (6,10) 

F26 = Net IR flux divergence in the layer (2,6) 

A workable scheme for making these computations is presented 
below for various combinations of cloud cover [CL(1), CL(2)]. 
For example, in order to compute 

Bi o 

F i o * = B 10 - B / 0 e wc (U,C,T)dB (3-1) 
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where B k is the Stef an-Boltzmann blackbody flux at T R , one 

must have a good representation of the emissivity e as a 

wc 

function of both water vapor and carbon dioxide absorber 
masses (U and C, respectively) along the sounding which is 
the path of integration. For the integral evaluation of 
Eq. (3-1), Sasamori (1968) proposed temperature-independent 
emissivity formulations e wc as exemplified for the absorber 
masses in layer (8,10) as follows: 

e (8,10) = {.240 logj o [U(8,10)]+.622}+.07262{ (1. 

W L/ 

- . 62556 [ U( 8 , 10 ) + . 0286] • 26 )[ logj „C(8,10) 

+1.064]} (3-2) 

Formulas of the type given by Eq. (3-2) apply for T >. 210°K. 
The quantity in the first brace of Eq. (3-2) indicates the 
water vapor emissivity alone and that in the second brace is 
the additional contribution by that part of the C0 2 emissiv- 
ity not already included in the water vapor emissivity ex- 
pression . 

Eq. (3-2) was used in the calculation of the atmospheric 
emissivities over all layers (6,10), (4 , 10) , . . . . , (1 , 10) , for 

which cumulative mass values U(6,10), C(6,10), etc., were 
available from the radiative soundings. Thus Eq. (3-2) 
proved very versatile. For completion of the integration 
from level k=l to the origin (B=0, e=1.0) of Yamamoto's Ra- 
diation Chart, the U(T ), C(T ) soundings were taken to be 

K K 

isothermal at T=Tj to the top of atmosphere (to level "0"). 
This leaves only the temperature-dependent emissivity e wc 
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[U(0,10),T 1 ] which applied for T <_ 210°K, also computed after 
Sasamori (1968) by means of 

l UCO.IO),!'! = 8 . 34T! ^ ‘ 353 lo 6i oU(0,10)-.44] 

V/ c 

X U(0,10) 1- * 03455 logioU(0,10)-.705] 

x [ 8 . 0/ ( . 353 logj 0 U(O,lO)+3. 56)] 

+ .07262{ (1.-. 62556 [U( 0,10)+. 0286] * 26 ) 

x [ logi 0 C(O,lO)+l. 064] } (3-3) 

The first term is temperature-dependent, as well as dependent 
upon the summed absorber masses relative to the reference 
level (in this example, level k=10). The term in the brace 
involves essentially the same form of the C0 2 correction to 
the total emissivity as seen in Eq. (3-2) except that it re- 
fers to the full layer (0,10). Eq. (3-3) is applicable with 
appropriate final parameters U(0,k), C(0,k), and T x to the 
determination of the downward flux contributions through 
level k for all of the reference level cases to be consid- 
ered below, i.e. k=2,6,10. 

B . NET FLUX F ] 0 * 

The radiative soundings were evaluated in the form of 
the parameters U(k,10), C(k,10), and T^, level by level. 

Also listed for each radiative sounding are the cloud cover 
parameters CL(1) and CL(2), which in general are non-zero. 

The grid area may be visualized as composed of areal frac- 
tions or weights defined as follows: 
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W(0,0) = ( 1-CL( 1 ) ) ( 1-CL( 2 ) ) : totally clear fraction 

W(1,0) = (CL( 1 ) ) ( 1-CL( 2 ) ) : overcast in upper layer only 

(3-4) 

W(0,1) = ( 1-CL( 1) ) (CL( 2 ) ) : overcast in lower layer only 
W(l,l) = (CL( 1 ) ) (CL( 2 ) ) : overcast in both layers 

The composite radiative effect of the true cloud cover 
over the grid area is then reconstructed by multiplying the 
results of the respective cases defined above by the appro- 
priate weight factors. The four cases will now be defined 
in order. 

1 . Net Flux F )n * with Clear Skies 

Defining the Stef an-Boltzmann blackbody flux at T^ 
as B k> where B^ = 1.170403 x 10“ 7 T^ ly(day) -1 , the net flux 
at the surface under the clear sky case is defined as fol- 
lows : 

F 10 *(O,O) = B 10 -.5{e wc (8,lO)(B 10 -B8)+(B 8 -B 6 )[e wc (6,lO) 

+ ewc (8,10)]+[ ewc (4 J 10)+e wc (6,10)](Be-B 4 ) 

(3-5) 

+(B 4 -B 2 )[e wc (2,10)+e wc (4,10)]+[e wc (l,10) 
+ ewc (2,10)](B 2 -B, )+e wc [(0,10),T 1 ]B 1 } 

2 . Net Flux F 10 * with Overcast Clouds in Upper Layer 
Only (k=4 to k-6 ) 

Fj o * ( 1 , 0 ) = (B 10 -B 6 ){l.-.5[e wc (8,lO)(B 10 -B 8 ) 

(3-6) 

+ (e wc (8,10) + £ : wc (6 ) 10))(B 8 -B 6 )]} 
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3 . Net Flux F in * with Overcast in Lower Layer Only 



(k=8 to k=9) 

Fio*(0,l) = (B 10 -B 9 )[l.-.5e wc (9,10)] (3-7) 

4 . Net Flux F 10 * with Overcast Clouds in Both Layers 
(4,6) and (8,9) 

F, 0 *(l,l) = (B 10 -B 9 )[l.-.5e wc (9,10)] (3-8) 

It should be noted that Eq. (3-7) and Eq. (3-8) are identi- 
cal . 

5. Composite Fi 0 * Calculations 

The composite results are obtained by multiplying 
the above defined four cases by their respective weights as 
defined by Eq. (3-4) to obtain 

F 10 * = [l-CL(2)]{(B 10 -B 6 )-.5e wc (8,lO)(B 10 -B 8 ) 

+[e wc (8,10)+e wc (6,10)](B 8 -B 6 )} 

+(l-CL(2))(l-CL(l)){B 6 -.5[(e wc (6,10) 

+e wc (4,10))(B 6 -B 4 )+(e wc (4,10)+e wc (2,10))(B 4 -B 2 ) 

+(e wc (2,10)+e wc (l,0))(B 2 -B 1 )+e wc ((0,10),T 1 )B, ] } 

+CL(2){(B 10 -B 9 )[l-.5e wc (9,lO)]} (3-9) 

It is clear that Eq. (3-9) reduces to F 10 *(O,O) when CL(1) = 
CL( 2 ) = 0, and to F 10 *(l,l) when CL(1) = CL(2) = 1.0. In a 
similar manner the other two cloud cases may be recovered by 
setting the weights to CL(1) = 1.0, CL(2) = 0 and CL(2) = 1.0, 
CL(1) = 0, respectively. 
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The above described procedure for obtaining the net 
flux can similarly be extended to levels k=6 and k=2 as 
listed in the following subsections. 

The composite F 10 * values will be compared in Section 
III.F with climatology (Malkus, 1962) to ensure the reason- 
ableness of this calculation. 

C. NET FLUX F 6 * 

The individual cases for F 6 * can be obtained from the 
composite F 6 * as described above for F 10 *. The numerical 
from using the Sasamori type emissivities [Eqs. (3-2) and 
(3-3)] for the composite F 6 * is as follows: 

F 6 * = [l-CL(l)]{B 8 -.5[e wc (6,8)(B 8 -B 6 )+e wc (4,6)(B 6 -B 4 ) 

+(e wc (4,6)+e wc (2,6))(B 4 -B 2 )+(c wc (2,6)+e wc (l,6))(B 2 -B 1 ) 

+e wc ((0,6),T, )B 1 ]}+(1-CL(1))(1-CL(2)){(B 10 -B 8 )[1. 

"•5( e wc (6.8)+e wc (6,10))]}+CL(l){(B 8 -B 6 )[l-.5(e wc (6,8))]} 

+CL(l)(l-CL(2)){(B 10 -B 8 )[l-.5(e (6,8)+e (6,10))]} 

WC WC (3-10) 

D. NET FLUX F 2 * 

In a similar fashion to F 10 * and F 6 * the net flux at 
level k=2 can be calculated from the radiation chart using 
the trapezoidal rule as described by Jenks (1974). Again 
the numerical formulation of the individual cloud cases can 
be obtained from the composite net flux formula as described 
in Subsection III.B.5. The composite F 2 * formula is 
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F 2 * = [l-CL(l)]{Ba-.5[e wc (2,4)(B 4 -B 2 )+(e wc (2,4) 
+ e wc ( 2, 6 ))(B6-B, )+(e wc (2,6)+e wc (2,8))(B 0 -B 6 ) 

( 3 — 1 1 ) 

+e wc (l,2)(B 2 -B 1 )+e wc ((0,2),T 1 )B, ] } +( 1-CL( 1 ) ) ( 1-CL( 2 ) ) 
x{(B 10 -B e )[l-.5(e wc (2,8)+e wc (2,lO)] }+CL(l){B„ 

■• 5 t £ wc (2>4)(B ^" B2)+e wc (1 ’ 2)(B2_Bl ) +g wc ((0 ' 2) ’ T i )B ^ } 

If the terms ^ wc ( (0, 2) ,T X )B X and e wc ( 1 > 2 ) (B 2 -Bj ) are set 
to zero in Eq. (3-11) a new parameter can be defined as the 
net outgoing long-wave radiation and labeled F 2 . This pa- 
rameter was computed to compare with satellite climatolog- 
ical observations over the oceans (Raschke, et al., 1973). 

E. APPLICATIONS TO HEAT BALANCE COMPUTATIONS 

For computation of the radiative "balance" at the top of 
the earth-tropospheric system only F 2 * of Eq. (3-11) and the 
total insolation absorbed below k=2 are required. At the 
surface, a radiative balance is computed at each gridpoint 
using the absorbed insolation and the net flux F 10 * of Eq. 
(3-9). Also at the surface a "convective" balance is com- 
puted by including the contributions of evaporation and sen- 
sible heat transport (see Sec. V). Net flux Fi 0 *, evaporation 
and sensible heat transport are all considered negative con- 
tributions at the surface. The atmospheric radiative bal- 
ance requires knowledge of the long-wave cooling effects 
caused by the flux divergences F26 and F610 as defined by: 

F26 = F 2 * - F 6 * ( 3-12a ) 

F610 = F 6 * - F ! o * ( 3-12b) 
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As will be seen in a later section, these parameters act as 
destabilizing influences in their respective layers when as- 
sessed in relationship to other processes in cross sectional 
display form in Section VI. 

F. COMPARISONS WITH PUBLISHED RESULTS 

Table II gives a visual comparison of the outgoing long- 
wave radiation, F 2 (as contrasted with F 2 *), found by this 
study as compared to the mean values published by Raschke 
et al. (1973). Also displayed are this study's global net 
flux, F 10 *, computed using Eq. (3-9) as compared to the an- 
nual mean deduced from Malkus (1962, Figs. 9,10). The 16 
October 1973 results have been interpolated, gridpoint by 
gridpoint, to 5° latitude increments for each of the four 
lines (see Fig. 2) and then latitudinally averaged to arrive 
at the "global" value for each latitude. The units of the 
values from Malkus (actually due to Budyko (1956)) have been 
converted from kcal cm -2 yr -1 to those of cal cm- 2 min -1 used 
in this study. The global mean values listed at the bottom 
of each column were derived by cosine weighting the latitu- 
dinal values listed in Table II (see also Sec. VII). 

It should be noted that the values of F 2 and Fi 0 * found 
by this study for latitudes 20°S, 10°S, and 60°N are averages 
of only two "sounding" lines while the other latitude values 
are averaged over all four meridional lines, the difference 
being due to the number of soundings available near the ex- 
tremities of the FNWC cross-sections. 
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The consistently smaller values of net flux, Fi 0 *, found 
by this study might be explained by the fact that Budyko ' s 
results are yearly climatological values, whereas the Fi 0 * 
values of this study, involving data from 16 October 1973, 
occurs near the time of the oceanic maximum of temperature 
and higher than annual mean cloudiness. This would reduce 
the net flux Fio* computed here. A latitudinal maximum net 
flux, Fio*, was observed around 40°N both for this study and 
the comparison figures of Malkus. At this location, Figs. 
7(b), 9(b), 10(b) show that there was a minimum opaque cloud 
cover for the time of this study. 

The comparison of outgoing long-wave fluxes, F 2 (Table 
11(b)), appears very consistent when comparing only northern 
hemispheric latitudes, considering that Raschke's values are 
based upon the mean for a month rather than that for a single 
day. Here maximum values of very nearly the same order oc- 
cur in the subtropics where minimum cloudiness occurs to im- 
pede the emission of outgoing terrestrial flux from the earth 
and lower atmosphere. There is in both this study and in 
Malkus', a well defined poleward gradient of F 2 by latitudes. 

The comparison of the computations of Fio* and F 2 made 
here with climatology (Table IIa,b) cannot confirm or deny 
the specific cloud-parameterization formulas (Eq. 2-13a,b) 
used in this study, but they indicate that they have a gen- 
eral reliability. 
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(a) 




(b) 




LAT 


F i o* 


Fi o * 


f 2 


f 2 




(Malkus ) 


(This Study) 


(This Study) (Raschke et 


20° S 


. 100 


. 068 


.333 


.40 


10°S 


. 090 


. 054 


.341 


. 39 


0° 


. 080 


.068 


.355 


. 36 


10°N 


. 090 


.079 


.378 


.34 


20°N 


.114 


.078 


.389 


.40 


30°N 


. 110 


. 081 


.388 


.37 


o 

o 


. 120 


. 107 


.376 


.34 


50°N 


.110 


. 086 


.346 


. 31 


60°N 


.120 


.067 


.315 


.29 


Global 

Mean 


.1018 


. 0757 


.3601 


. 3615 



Table 11(a). Comparison of net flux at the surface, F 10 *, 
as found by this study for 16 October 1973 and the annual 
mean F 10 * derivable from Malkus (1962, Figs. 9,10). 

(b). Comparison of net flux to space, F 2 , as found 
by this study and by Raschke et al. (1973), based upon 
NIMBUS III heat budget studies. 
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IV. SOLAR RADIATION 



A. COMPOSITION OF SOLAR INSOLATION 

At the top of the atmosphere (k=0) this study assumed a 
solar constant of 2.00 ly min -1 (Joseph, 1971); furthermore, 
this constant was assumed subject to a four percent attenua- 
tion above the tropopause due to ozone and oxygen. Thus the 
effective solar constant at level k=2 in this study is 1.92 
ly min -1 . 

To compute the effective solar insolation at the tropo- 
pause the following formula was used 

F(2) = S[j-]~ 2 Cos Z (4-1) 

m 

where S = effective solar constant at level k=2 

Cos Z = cosine of the zenith angle, a function of 
the Julian date 

r/r = ratio of the actual earth-sun distance to 
the mean earth-sun distance, a function of 
the Julian date. 

The Smithsonian Meteorological Tables (List, 1958) list the 

ratio r/r for 15 October at 0000GMT as 0.997165, this value 
m 

was utilized for all four data lines. The cosine of the 
zenith angle was determined by 

Cos Z = Sin 4> Sin 6 + Cos <t> Cos 6 Cos h (4-2) 
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where 



<P = latitude 

6 = solar declination for 15 October = -8° 13.5' 

(Table 169, List, 1958) 

h = hour angle 

The value of sine <j> was calculated using one of two differ- 
ent formulae, depending on the data-line used for the compu- 
tations, in terms of the FNWC map coordinates (I,J) as in 
Eq. (4-3a,b), conversely for these lines one may solve for I 
in terms of sine <}> as in (4-3c,d): 



Lines 1, 


3,4 


Sin <j) = 


973. 

973. 


.752 

.752 


- 2(32-1 ) z 
+ 2 ( 32- 1 ) 2 


(4- 


-3a) 


Line 2 




Sin <j> = 


973. 

973. 


.752 

.752 


- ( 32-1 ) 2 
+ ( 32- 1 ) 2 


(4- 


-3b) 


Lines 1, 


3,4 


I = 32 


- 22. 


.065 


r Cos <f> , 

l l + Sin J 


(4- 


-3c) 


Line 2 




I = 32 


- 31. 


.205 


, Cos <(> , 

l l + Sin <f> J 


(4- 


-3d) 






. . ,25 




for 


Lines 1,2 






I = 


< 8,. 


. . ,25 




for 


Line 3 








( 63, 


. . . , 38 




for 


Line 4 







Here I is the abscissa distance on the FNWC grid, Fig. 2, 
and varies by line as described in Sec. II. The soundings 
for' lines 1, 2, and 3 were all taken at OOOOGMT with the so- 
lar noon existing at the 180 meridian, therefore the hour 
angles for these three lines were 55° , 10° , and 35° , re- 
spectively. For line 4, the soundings were taken 12 hours 
earlier with solar noon at the Greenwich meridian, giving an 
hour angle for line 4 of 35°. 
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The composition of solar insolation was partitioned in 
this study, in accordance with Joseph (1971), by dividing 
the insolation, F(2) at level 2, into two parts, one of 
which was subject to water vapor absorption but not to Ray- 
leigh scattering, F(A), and one which was only subject to 
Rayleigh scattering, F(S), in a clear moist atmosphere. 

These partitions may be easily computed as follows: 

F(A) = 0.349 F(2) (4-4) 

F(S) = 0.651 F(2) (4-5) 

Cloud-reflectivity is inherent in both F(A) and F(S) solar 
energy portions with the presence of two cloud decks. 

To further define the solar insolation partitions, the 
F(A) energy consists of those wavelengths X > 0.9 ym where 
absorption by water vapor and carbon dioxide bands is the 
dominant attenuation process for clear air. The F(S) energy 
is composed of all wavelengths X _< 0.9 ym where water vapor 
absorption is considered negligible. 

B. DISPOSITION OF F(S) INSOLATION 

As described by Jenks (1974), there are two reflective- 
type parameters which must be considered when studying the 
energy contained in the shorter wavelengths, F(S). The 
first parameter is the Rayleigh clear-sky albedo, (after 
Coulson, 1959), which was empirically formulated using the 
least squares best-fit technique by Joseph (1971) to be 
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(4-6) 



a(E) = .085 + .25074 [ logx 0 Sec Z)] 

r o 

P 0 = 1013.25 mb 

Since mean sea level is very nearly 1000 mb, this study as- 
sumed that 7 t/P 0 -1. The value of Sec. Z is the inverse of 
Eq. (4-2). 

The second reflective parameter considered is the surface 
albedo, which is a function of the cosine at the zenith 
angle over ocean surfaces. The surface albedo, as utilized 
by this study, follows Gates et al. (1971). 

a (G) = max (.06, .06 + .54 (.7 - Cos Z)} (4-7) 

As described in Sec. Ill, four distinct cloud cases are pos- 
sible with a two-layer cloud model. The disposition of F(S) 
under each of these cases will be discussed in the remainder 
of this subsection. 

1 . Clear Sky Case 

With clear sky the F(S) insolation is modified only 
by Rayleigh sky reflectivity, a(R), and surface reflectivity, 
a(G), as given by Eqs. (4-6) and (4-7). Jenks (1974)points 
out that with multiple reflections possible between the 
earth and atmosphere, the clear sky insolation absorbed at 
the earth's surface (k=10) after scattering is given by 

IS10(0,0) = F(S) [l-a(R)] [l-a(G)]/[l-a(R)a(G)] (4-8) 

where the notation (0,0) denotes CL(1) = CL(2) = 0.0 as dis- 
cussed in section III. 



43 



2 . Cloudy Sky Cases 



With the presence of cloud layers, an additional in- 
solation modifier must be considered, cloud reflectivity. 

In this study the values suggested by C.D. Rodgers (1967) 
were adopted, namely, R(l) = 0.54 for the higher clouds (be- 
tween k=4 and k=6) and R(2) = 0.66 for the lower clouds (be- 
tween k=8 and k=9) . 

Consider first the case (1,0) where the fraction of 
the grid area has overcast clouds only at the high level, 
that is where CL(1) = 1.0, CL(2) = 0.0, and R(2) = 0.0 so 
that the earth-absorbed F(S) insolation is 

IS10( 1 , 0 ) = F(S)(1-R(1))( 1-a ( G) ) / ( 1-R( 1 )a ( G ) ) (4-9) 

The F(S) insolation absorbed by the surface above 
which there is present only a lower cloud overcast, [denoted 
by (0,1)] where CL(1) = 0.0, CL(2) = 1.0, and R(l) = 0.0 is 
given by 

IS10( 0,1) = F(S)(1-R(2))( 1-a (G) )/ (l-R(2)a(G) ) (4-10) 

If overcast clouds are present in both layers, then 
the conditions CL(1) = CL(2) = 1.0 [or in abbreviated symbol- 
ism, case (1,1)] applied, and reflectance also exists between 
the cloud layers as well as between clouds and space. Thus 
the absorbed F(S) insolation at the earth's surface as for- 
mulated after Arakawa (1972) is given by 

IS10(1,1) = F(S)[(l-R(l))(l-R(2))(l-a(G))]/[l-((R(l)R(2) 
+R( 2)a(G)+R(l)a(G) )+2(R(l)R(2)a(G) )] (4-11) 
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a formula which resembles very closely a combination of Eqs. 
(4-9) and (4-10). 

Thus far in this subsection the sky was considered 
either overcast or clear at each level without regard for 
the actual cloud conditions CL(1) and CL(2). 

3 . Composite F(S) Insolation 

If the areal coverage weights, as defined by Eq. 
(3-4), are applied to the appropriate right sides of Eqs. 
(4-8), (4-9), (4-10), and (4-11), then the composite surface- 



absorbed insolation, IS10, results 

IS10 = IS10(0,0)W(0,0)+IS10(1,0)W(1,0) 
+ IS10(0,1)W(0,1)+IS10(1,1)W(1,1) 



(4-12) 



It is apparent that F(S) insolation reflected to space may 
be found by subtracting IS10 from F(S). 



C. DISPOSITION OF F(A) INSOLATION 

This subsection deals with that portion of the solar in- 
solation which is absorbed by atmospheric water vapor and 
carbon dioxide, namely the longer wavelengths. The discus- 
sion will again be divided into the four different cloud- 
layer cases. 

1 . The Clear Sky Case 

Formulation enabling the determination of absorption 
values for particular layers of interest was provided by use 
of the Manabe-Moller absorptivity function which may be gen- 
eralized as follows 

a( 2 , 2k ) = ■. 271[U(2,2k) Sec Z]* 3 0 3 ( 4-13 ) 
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Here absorptivity a is applied over the pressure-scaled water 
vapor mass between levels 2 and 2k (Fig. 1) along the slant 
path of the zenith angle Z. Recalling that the water vapor 
mass above level 2 was assumed negligible, the absorbed solar 
insolation, for example, in the layer (2,6) is written as 

A( 2 , 6 ) = 0.271 F(A)[ U(2 , 6) Sec Z] ,3 ° 3 (4-14) 

and the absorbed solar insolation in the layer (6,10) is 
given by 

A(6, 10) = A( 2 , 10 ) - A(2 , 6) (4-15) 

For a schematic depiction of the F(A) insolation disposition 
in the clear sky case [see Fig. 3 in Jenks (1974)] . 

To obtain that portion of F(A) insolation which is 
absorbed at the earth's surface, the atmospheric transmissiv- 
ity (l-a(2,10)) and the surface transmissivity after reflec- 
tance are combined (by the product-law of transmissivities) 
to give IA10 as follows: 

IA10( 0 , 0) = F(A){1-.271[U(2,10) Sec Z] * 3 0 3 } ( l-a(G) ) 

(4-16) 

Here a(G) is defined by Eq. (4-7) 

In the remainder of this subsection, which discusses 
the cloudy layer cases, representative cloud reflectivities 
and cloud absorptivit ies were adopted, after C.D. Rodgers 
(1967), for the two possible cloud layers. The cloud re- 
flectivities were RA(1) = .46 and RA(2) = .50 while the 
cloud absorptivities were A(l) = .20 and A(2) = .30. Note 
that the reflectivities for the F(A) wavelengths differ from 
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those adopted for the F(S) wavelengths. The procedure of 
considering the cloud conditions to be overcast whenever 
they appear and then applying the appropriate weighting fac- 
tors in the composite summation will again be followed as in 
Sec. IV. B. 

To simplify the discussion for the cloud-covered 
cases, the (1,1) case will be presented first as it contains 
representative type equations for the remaining two cases, 
(1,0) and (0,1). 

2 . Overcast in Both High and Low Cloud Layers 

The following set of formulas detail the disposition 
of incoming solar insolation (F(A)) from level k=2 to the 
earth's surface to permit determination of the amount of in- 
solation absorbed by the atmospheric layers and by the earth's 
surface. The dashed separation lines are introduced to sub- 
divide the absorption and reflection physics of the model 
into subsections which permit the analysis to proceed more 
or less within successive 200 mb layers (see also Fig. 3). 

F4i = F(A)(1-.271[U(2,4) Sec Z]' 303 ) 

F4+ = F41 (RA( 1) ) 

(4-17a) 

A241 = F(A) . 271 [U(2 , 4) Sec Z]- 303 

A24t = F4t . 271 [ U( 2 , 4 ) Sec Z]* 303 

F61 = F4+(1-RA(1)-A(1)) 

A46 = F41 ( A( 1 ) ) (4-17b) 

A26 = A244- + A24t + A46 
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Figure 3. Schematic representation of F(A) insolation 
disposition in the case of two overcast layers. 



48 



TD68 = 
TD910 


= 1-. 271[U(6, 8)5/3] * 303 

(4-18) 

= 1-. 271[U(9, 10)5/3] * 303 



F84- = F6t (TD68 ) 
F8t = F84- (RA(2 ) ) 
F6t = F8+ (TD68 ) 



F6+ i = 


= F6t ( RA( 1 ) ) 


F84-4- = 


= F6tt (TD68 ) (4-19a) 


A684- = 


= F6+ - F8^ 


A68t = 


= F8t — F6+ 


A684'4 i 


— FG'V 4 r — F84 , 4' 


A68 = 


ASS'I' A68+ + A68i + 


F9t = 
A89 = 


F8t ( 1-RA( 2 ) - A(2) ) + F8tt(l-A(2)) 

( 4-19b ) 

F8t(A(2)) + F8tt ( A(2 ) ) 


FlOt = 


= F9t (TD910) 


FlOt = 


= FlOt (a (G) ) 


F9t = 


FlOt (TD910) 


F9t t : 


= F9t (RA(2) ) 


FlOt 1 


= F9t 4- (TD910) (4-19c) 



A9104- = F94- - F104- 
A910+ = FlOt - F9t 
A9104- 4- = F94-4- - FlO+t 



A910 : 


= A910t + A910t + A910t t 


A810 : 


= A89 + A910 
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II 










In the preceding equations the symbol "A" followed 
by a set of digits refers to the F(A) insolation absorbed in 
the indicated layer, for example, A26 indicates the absorp- 
tion in the layer k=2 to k=6. The symbol "F" followed by a 
k value indicates the amount of insolation passing through 
the particular k level. The vertical arrow modifiers indi- 
cate the direction of insolation passage, i.e., 1 downward 
insolation, t upward reflected insolation, and -j-i downward 
reflected insolation. Terms involving the symbols "TD" as 
expressed by the functions of Eq. (4-18), indicate the 
Manabe-Moller transmissivities for insolation beneath a 
cloud layer. In the latter situation, the term Sec Z in 
Eq. (4-13) is effectively replaced by the mean slant-path 
Sec Z ^ 5/3 (Katayama, 1966). 

Note that the effect of multiple reflections between 
clouds or between the earth's surface and the lower cloud 
has been incorporated to the extent of two reflections, with 
the lowermost reflecting surface absorbing the remaining im- 
pinging insolation. Computations indicated that the insola- 
tion remaining after two reflections was too small to warrant 
further reflections. Also insolation reflected upward from 
a lower interface (cloud or ground) to the base of an upper 
cloud deck has not been subjected to cloud absorption. This 
tends to reduce very slightly the secondary cloud-absorption. 

Note finally that between any two reflecting surfaces 
there are three contributions to absorption of solar insola- 
tion within the layer. These contributions are denoted by 
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the arrow-identifiers which indicate the segment of path 
being traversed. 

From Eqs. (4-19c) the impinging F(A) insolation at 
the earth's surface can be expressed 

TRANA( 1,1) = FlOi + F101I (4-20) 

The F(A) insolation which is actually absorbed by the earth's 
surface (see Fig. 3) may be written as 

IA10( 1,1) = F101 ( 1-a (G) ) + F104- 1 (4-21) 

3 . Disposition of F(A) Insolation with an Upper Overcast 
With a cloud layer present only in the upper layer 
the equations depicting the disposition of incoming insola- 
tion becomes a simplified subset of the previous case: 

A241 = F(A) . 271[U(2 , 4) Sec Z]' 303 

A24t = F(A){1-.271[U(2,4) Sec Z]’ 303 }RA(1) (4-22a) 

x{ . 271 [ U( 2 , 4 ) Sec Z] * 3 0 3 } 

A24 = A24+ + A24t 

F6I = F(A){1-.271[U(2,4) Sec Z] * 3 0 3 } [ 1-RA( 1 )-A( 1 ) ] 

A46 = F(A){1-.271[U(2,4) Sec Z]- 303 }A(1) (4-22b) 

A26 = A24 + A46 
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F104- = F6-KTD610) 

FlOi = FlOKa(G)) 

F6 f = F10t(TD610) 

F64-4- = F6 1 (RA( 1 ) ) 

F104-4- = F64--KTD610) (4-22c) 

A6104- = F61 - FlOf 
A610t = FlOt — F6+ 

A610-H- = F6 + 4- - FIO-H- 

A610 = A610+ + A6104- + A610 + 4- 

The variables used above are defined in a similar manner to 
those in the (1,1) case. The impinging insolation at the 
earth's surface may be formulated as follows 

TRANA( 1 , 0) = F10+ + F10+1 (4-23) 

while the F(A) insolation absorbed by the earth in this case 
(Fig. 4) is given by 

IA10( 1,0) = F101 ( l-a(G) ) + FlOii (4-24) 

4 . Disposition of F(A) Insolation with a Low Overcast 
With an overcast lower cloud layer the disposition 
formulas are as follows: 
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Figure 4. Schematic representation at F(A) insolation dis- 
position with an upper overcast layer only. 
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A261 = 


= F(A) .271 (U(2 , 6) Sec Z] • 3 0 3 


A684- = 


= F( A) . 271{ [ U( 2 , 8) Sec Z] * 3 0 3 -[ U( 2 , 6 ) Sec Z] • 3 0 3 } 


F8| — 


F(A){ 1-.271[ U(2,8) Sec Z] • 3 0 3 } 


F8t = 


F8t [ RA(2 )] 


F6t — 


F8t{l-.271[U(6,8) Sec Z] * 3 0 3 } (4-25a) 


F2+ = 


F8t{l-.271[U(2,8) Sec Z]- 303 } 


A68t = 


= F8t — F6t 


A261 = 


= F6t - F2t 


A26 = 


A26I + A261 


F9+ = 
A89 = 


F8+[ 1-RA(2) - A(2)] 

(4-25b) 

F8l [ A ( 2 )] 


F104- = 


= F9| (TD910) 


FlOt = 


= F104- (a (G) ) 


F9t = 


F10+ (TD910) 


F94-4- = 


= F9t (RA( 2 ) ) 


F104-4- 


= F94-4- (TD910) (4-25c) 


A910+ 


= F94- — F104- 


A910+ 


= FI Of - F9t 



A910++ = F91 4- - F10+ + 



A910 = 


= A910+ + A910t + A910I4- 


A610 = 


= A684- + A68t + A89 + A910 



Again the variables are defined as stated before in the 
(1,1) case. Likewise the incident flux is defined at the 
earth's surface as 

TRANA(0, 1) = F10 + + F10++ (4-26) 
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while that portion which is absorbed by the earth, (Fig. 5) 

is 

IA10( 0 , 1 ) = FlOt ( 1-a (G) ) + F10W (4-27) 

5. Weighted F(A) Layer-Absorptions and Surface-Absorp- 
tion Insolation 

As described in Section IV. B. 3 on the composite F(S) 
insolation, the values of the absorbed F(A) insolation in 
the layers and of the earth's surface which were calculated 
for the four cloud cases (0,0), (1,0), (0,1), and (1,1) must 

be appropriately weighted by the values from Eq. (3-4) to 
obtain the composite F(A) layer- and surf ace- absorbed inso- 
lation . 



A26 = A26 (0>0) W(0,0)+A26 (1>0) W(1,0)+A26 (0>1) W(0,1) 



(4-28) 



+A26 (11) W(1,1) 



A610 = A610 (0 0) W(0,0)+A610 (1 q) W( 1 , 0)+A610 ( Q , ± } W( 0 , 1 ) 



(4-29) 



+A610 (1 1} W(1,1) 



IA10 = IA10(0,0)W(0,0)+IA10(1,0)W(1,0)+IA10(0,1)W(0,1) 

(4-30) 

+IA10(1,1)W(1,1) 



6. Absorptivity (ABA) by Layers 

It was convenient to work with the absorptivity as 
well as the actual insolation energies absorbed by each lay- 
er and at the surface. To obtain the absorptivity, the 
total incoming insolation at k=0 must be known. This value 
called FADJ is formulated as follows 
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Figure 5. Schematic representation of F(A) insolation 
disposition with a lower overcast layer only. 
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(4-31) 



FADJ =2.00 (r/r ) -2 Cos Z 
' m 

By taking the ratio of the insolation absorbed in the tropo- 
sphere to that incident at the top of the atmosphere, the 
full tropospheric absorptivity is formed 



ABA 



A26 + A610 
FADJ 



(4-32) 



Similar formulas may be expressed for the partial absorptiv- 
ities for the layers (2,6) and (6,10). 

D. REFLECTION OF F( A) AND F(S) ENERGY 

To determine the albedo and fractional absorptivity of 
the earth's surface, both the F(A) and F(S) fractions of so- 
lar insolation must be considered. 

1. Albedo 

The amount of F(A) insolation which is reflected 
back to space is determined by the following 



REFA = F( A) - A26 - A610 - IA10 . (4-33) 

here A26, A610, and IA10 are defined by Eqs. (4-28), (4-29), 
and (4-30), respectively. 

Because the F(S) wavelengths do not undergo atmos- 
pheric absorption, the reflected F(S) insolation was easier 
to calculate 



REFS = F(S) - IS10 (4-34) 

IS10 was defined by Eq. (4-12). By combining Eqs. (4-33) 
and (4-34) the total reflected insolation to space is de- 
fined by 
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REF = REFA + REFS 



(4-35) 



The resulting REF values were calculated at each gridpoint 
and are meridionally displayed by latitudinal interpolation 
in Figs. 7(a) 10(b) . 

By taking the ratio of the total reflected energy to 
the total incoming energy the earth-atmospheric system albedo 
is defined 

ALB = REF/FADJ (4-36) 

where FADJ is defined by Eq. (4-31). 

2 . Composite Absorptivity (ABG) by the Earth's Surface 

The fractional absorptivity of the earth's surface 
was determined by the ratio of the total insolation absorbed 
by the earth to the total incoming insolation. The total 
insolation absorbed by the surface is determined by the prop- 
erly weighted sums of IA10 and IS10, defined by Eqs. (4-30) 
and (4-12), respectively. 

ABG = ( I A10 + IS10) /FADJ (4-37) 

The result of Eq. (4-37) gives the weighted absorptivity of 
earth. It is also of interest to determine the total inso- 
lation incident at the earth's surface just before absorp- 
tion by earth. Th'e result is given by the properly weighted 
sum of terms of the form 

TRAN = TRANA + [ IS10/ ( 1-a (G) ) ] (4-38) 
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It should be noted that each of the four terms of form IS10 
(Eq. 4-12) had a factor (l-a(G)) in the numerator so that the 
transmitted F(S) insolation incident at earth before absorp- 
tion need only be divided by (l-a(G)). The parameter TRAN 
is a simulation of the insolational energy incident at a 
pyrheliometer . Note finally that the ratio 

ATRAN = TRAN/FADJ (4-39) 

represents the transmissivity of the troposphere. 

3 . Computational Check 

Recalling the assumption of a four percent attenua- 
tion of the incoming solar insolation in the stratosphere, 
the summation of the fractional absorpt ivit ies ALB, ABA, and 
ABG should equal 0.96. This value provides a good computa- 
tional check for this model at each gridpoint and was in 
fact satisfied over the entire grid. 

E. STATISTICAL TESTS 

To further verify the foregoing calculations, statisti- 
cal computations were performed on several key parameters 
utilizing linear regression techniques from the BMD set of 
statistical computer programs (Dixon, 1973). 

1 . Clear Sky Case Regressions 

Three linear regressions were computed for the clear 
sky case, the first of which was the clear sky albedo, 

ALB( 0 , 0 , Z ) . The best-fit equation, using the secant of the 
zenith angle, Sec Z, as the predictor, was computed to be 

ALB( 0 , 0 , Z ) = 0.00650 + 0.09052 Sec Z (4-40) 
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with a regression coefficient R = .964. This highly signif- 
icant correlation coefficient confirms the model results 
which indicate that Sec Z is a good predictor of clear-sky 
albedo. In this connection, it should be recalled that a(G) 
and a (Ft) were modeled to be functionally related to Sec Z. 
The mean global clear-sky albedo was 

ALB( 0,0) = 0.164 

which corresponds to the mean secant of the zenith angle, 

Sec Z = 1.743, for this study. The statistics show that the 
largest contributions to error in the fit by Eq. (4-40) oc- 
curred at largest zenith angles, which occur at high lati- 
tudes . 

The two other clear-sky linear regressions performed 
were ABA(0,0,M) and ATRAN(0 , 0 , M) , where the predictor, M, 
was the water-vapor mass path length given by 

M = (U Sec Z)^[ log! o (U Sec Z) 4 ] (4-41) 

This parameterization of the water-vapor mass is formulated 
after Hanson (1967, 1971), who used a similar absorber-mass 
predictor in his empirical equations for ABA with both clear 
and partly cloudy skies. The resulting clear-sky regression 
equations in this study for ABA and ATRAN were 

ABA ( 0 , 0 , M ) = 0.10320 + 0.05535M (4-42) 

with R = .9925 and 

ATRAN ( 0 , 0 , M ) = 0.77494 - 0.04823M (4-43) 
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with R = .7142. The sample means of the statistics which 
entered into (4-42) and (4-43) are 

ABA ( 0 , 0 , M ) = 0.146 

ATRAN ( 0 , 0 , M ) = 0.738 

M = 0.776 (gm cm" 2 )* 

based on the original 93 cases. 

2 . Cloudy-Sky Case Regressions 

In addition to the clear-sky fractional amounts of 
ALB, ABA, and ATRAN, linear regressions were formulated for 
the composite cloudy-sky case, i.e. ALB[CL(1) ,CL(2)] , 

ABA [CL ( 1 ) , CL( 2 ) ] , and ATRAN [ CL( 1 ) , CL( 2 ) ] , respectively. The 
regression relationships between the clear- and cloudy-sky 
cases were determined for the same zenith angle at gridpoint 
(I,J) using the total opaque cloud cover (CL) as the predic- 
tor. The total opaque cloud cover approximates the effective 
cloud cover existing in a two-layer cloud model and is ex- 
pressed as 

CL = CL( 1 ) + CL( 2 ) - CL( 1 )CL(2 ) (4-44) 

The relationship between the cloudy-sky albedo and 
the corresponding clear-sky calculation was obtained by 
least squares with the following results 

ALB [CL ( 1 ) , CL( 2 ) ] = ALB( 0, 0) [ 1+2 . 09929CL+0. 84608CL 2 ] (4-45) 

and R = .9292. The means of the parameters used in Eq. (4-45) 
are 



61 



ALB[ CL( 1) ,CL(2)] = 0.406 



ALB( 0,0) = 0.164 
CL = 0.662 

The dependence of Eq. (4-45) upon CL is indicated by 
the high correlation coefficient (R=.9292). The global mean 
albedo for this study of 0.406 is somewhat higher than that 
reported in recent literature (ALB=.300) by Raschke et al. 
(1973). However, Raschke' s figure is an annual mean which 
includes the somewhat smaller cloud amounts over land. 

The tropospheric absorptivity of air in the composite 
cloudy case was developed analogously to (4-45) by least 
squares and gave the form 

ABA [ CL ( 1 ) , CL ( 2 ) ] = ABA( 0 , 0) [ 1+0 . 22305CL+0 . 17261CL 2 ] (4-46) 

with R = .927 and mean statistics for the above parameters 
ABA[ CL( 1 ) , CL( 2 ) ] = 0.180 
ABA( 0,0) = 0.146 
CL = 0.662 

The direct relationship between solar absorptivity and cloud 
cover is in general agreement with the results of Plante 
(1972) who found increased absorptivity with increasing CL. 
This relationship could be the reason for the recent reduced 
estimates of global albedo [Raschke et al. (197 o), Vondei 
Haar and Suomi (1971)]. The significant value of R indi- 
cates the degree of sensitivity between cloud amounts and 
atmospheric absorptivity. 
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The final statistical regression equation formulated 
deals with the cloud-sky transmissivity of the troposphere 
relative to the clear-sky calculation. Again the total 
opaque cloud cover has been used as the independent variable 

ATRAN [ CL( 1 ) , CL( 2 ) ] = ATRAN(0, 0) [1-0. 51477CL-0. 14646CL 2 ] 

(4-47) 

Here the correlation coefficient R = .997 and the mean sta- 
tistics over the 93 sample soundings are 

ATRAN [CL(1) ,CL(2) ] = 0.543 

ATRAN(0, 0) = 0.738 

CL = 0.662 

Note the high degree of significance, (R = .997) of Eq. 
(4-47), which formulation agrees with empirical results of 
Quinn (1971) and numerous other workers cited by Quinn in 
his 1971 paper. 

In general, the solar-model statistical relation- 
ships developed in this study show consistency with empiri- 
cally reported observations, which have already been noted, 
when consideration is taken for seasonal effects and the 
latitudinal limits of this study. 
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V. SENSIBLE AND LATENT HEAT TRANSPORT 



AT THE SEA-AIR INTERFACE 
A. GENERAL CONSIDERATIONS 

For a complete discussion of the heat balances of the 
earth and atmosphere for 16 October 1973, the vertical 
transports of sensible and of latent heats across the sea- 
air interface must be considered. For this purpose, the 
turbulent flux models for the planetary boundary layer, as 
used in the FNWC primitive equation model as described by 
Kesel and Winninghoff (1972) and Kaitala (1974) , was adapted 
here virtually without change. 



B. EVAPORATION 

The determination of the surface latent heat flux, E, 
over the ocean is based upon its large scale parameteriza- 
tion, principally after Langlois and Kwok (1969) 



E = Lp 



K* 



1 o 



e -e 

l+a* ( 77 ^ — 77 ^) 



q x " q 9 
Z 9 -Z ! 0 



(5-1) 



Z 9 — Z 1 0 



Here L is the latent heat of vaporization given by 
L = 596.73 - .601 T 10 (cal/gm) 

The other parameters in (5-1) are formulated as follows 



P 1 0 - P 1 0 / 1 



P — T ( 1000 \ 2/7 
° 9 i9( 900 ' 



e s = t 1 0 
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Z 9 -Zj o 




9 1 0 



T 



In 



1000 

900 



q x -q, = .62197 ^ -q(900) 



Note that at level k=10, the assumption P t 0 - 1000 mb is con- 
sidered valid and e is the vapor pressure at the top of the 
constant flux layer. Also q(900) is the mixing ratio at k=9 
as defined by Eq. (2-4). In Eq . (5-1) K* and a* are turbulent 
transfer coefficients having values 2 x 10 5 cm 2 sec -1 and 
5 x 10 1 * cm/°K, respectively, according to Kaitala (1974). 
However, values of the eddy coefficient, K* , of 2 x 10 5 and 
1 x 10 s were both tested in this study and it was decided 
primarily to use only the latter value. Furthermore, it is 
felt by this author that upon further testing, with a broader 
data base, it may be feasible to tune K* between the two 
values tested so as to make it dependent either upon latitude 
or season, or other pertinent variables. 

C. SENSIBLE HEAT TRANSPORT 

Sensible heating is computed as a function of the tem- 
perature difference between the surface, Tjo, and the air 
near the surface, T , which is at the top of the constant 

X 

flux layer. Since wind values were not carried in this 
study, the technique used by Mintz and Arakawa, as described 
by Kaitala (1974) and by Kesel and Winninghoff (1972), was 
first attempted for sensible heating 
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(5-2) 



Hr, — P 1 0 C ( 



K* 



°9-°S 

l +a+( _J__s } 



e -6 

){y„ - 9 x 



Z 9 — Z i 0 



} 



z 9 -z 



i o 



Here C - .239 cal/(gm°K), and, since the constant flux lay- 
er was assumed to be quite thin, G =T . v is a critical 
lapse rate (3 * 10 -5 °K/cm), and the remainder of the terms 
are defined as in the discussion following Eq. (5-1). The 
micrometeorological equivalent of Eq . (5-2) is given by 
Kaitala (1974) as 



H r = e-VWVV < 5 - 3 > 

It can be seen from Eqs. (5-2), (5-3) that the temperature 

lapse, T -T , for .the thin layer between the surface and the 
g x 

top of the constant flux layer becomes very critical in de- 
termining the amount of sensible heating (cooling) for each 
sounding location. 

In Eq. (5-2), the expression within the brace determined 
the sign of the sensible heat flux in the vertical. In ap- 
proximately one out of three cases, the braced expression 
gave rise to negative values of Hp , whereas the correspond- 
ing latent heat flux, E, of Eq . (5-1) was systematically 

positive. Based on Fig. 15 in Malkus (1962), which shows Hp>0, 
there seems to be no valid basis for having an erratically 
downward distribution of sensible heat over the open oceans. 
Consequently, this study assumed that E was properly as- 
sessed by Eq . (5-1), and that a reasonable set of values of 

Hp could be obtained by means of 
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I 






rE 



(5-4) 



Here the Bowen ratio r was determined from a set of best es- 
timates as a function of latitude. 

Using 20 values of r and their functionally related lati- 
tudes, taken from published results of Sellers (1965), 

Budyko (1956), and Jacobs (1951), a linear regression was 
best-fitted to | Sin <j> | to give 

r = 0.10910 - .41996 | Sin <t>\ + .89645 | Sin 4* | 2 (5-5) 

with a correlation coefficient of R = .982. Thus, using 
Eqs. (5-4) and (5-5), the sensible heat transport was ob- 
tained at each gridpoint, while always maintaining proper 
sign sense with E which was virtually always positive. 

D. COMBINED TURBULENT HEAT FLUX ACROSS THE SEA-AIR INTER- 
FACE 

To get a better feel for the effect that sensible heat- 
ing has at the surface, the summed terms E and Hp were uti- 
lized as heat-loss quantities in the total budget at the 
surface, therefore in terms of large scale "convective" dis- 
position 

E + Hp = E( 1+r ) (5-6) 

where r ranges in magnitude from approximately .05 to .5 
with increasing latitude. 
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VI . MERIDIONAL CROSS-SECTIONAL DEPICTION OF 
THE HEATING MODEL COMPUTATIONS 

A. GENERAL DESIGN OF HEAT-MODEL OUTPUT 

In this section, the computational subroutines already 
developed (Secs. Ill, IV, and V) are combined to perform a 
single time-step of the heating program developed for adapta- 
tion to the FNWC prediction model. The results represent 
computations of the relevant heating functions at each grid- 
point of the four meridians (Fig. 2). 

For the three Pacific cross-sections, the FNWC gridpoint- 
processed analyses for OOOOGMT, 16 October 1973, have been 
used, whereas for the 35°W (Atlantic) meridian, the FNWC 
analyses from 1200GMT, 15 October 1973, were utilized. For 
simplicity the heating model computations will be referred 
to as those pertaining to 16 October 1973. 

Fig. 6 depicts the key to the computations made at (I,J) 
gridpoint soundings, shown in the form of Table 1(b), and in 
effect represents the flow diagram of computations proceed- 
ing from top-to-bottom at each gridpoint. Figs. 7, 8, 9, 
and 10 show that the gridpoints considered are at 5° lati- 
tude increments, but this has been done merely for simplifi- 
cation of the interpretations and for comparisons with 
climatological data. In order to deduce a heating pai'ameter, 
Q( I ) , at a non-integral grid-value I corresponding to lati- 
tude 4> (in whole multiples of 5°), a Lagrangian cubic inter- 
polation scheme was utilized of the form 
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( 6 - 1 ) 



om = Qi ( i~2 ) ( 1-3 ) ( 1-4 ) + Q 2 (I-l)(I-3)(I-4) 

^ J ( 1-2 ) ( 1-3 ) ( 1-4 ) (2— 1)(2— 3)( 2-4 ) 

Q, (I-l)(I-2)(I-4) Q„ (I-l)(I-2)(I-3) 

(3— 1)(3— 2) ( 3-4 ) (4-l)(4-2)(4-3) 

Here I is determined by an equation of the form (4-3c), the 
points 1,2,3, and 4 are any four successive integral values 
of I with I lying in general between points 2 and 3, although 
the subroutine also works well if I is non-centered in the 
I-array . 

B. CROSS-SECTIONAL DATA AT LEVEL k=2 

The terms explained here are those occurring at the lev- 
el k=2. The discussion of all insolational parameters dis- 
cussed previously in section IV pertained to the specific 
time of day corresponding to the hour angle h at the instan- 
taneous times under consideration. Thus the incident inso- 
lation was given by Eq. (4-1), i.e. 

F(2 ) = S (-p-) -2 Cos Z 
m 

with Cos Z given by Eq. (4-2). For the specific map times 
the local solar times gave h=35° , 10° , 55° , and 35° for 
cross-sections 1,2,3, and 4 depicted in Figs. 7, 8, 9, and 
10 respectively. 

For each specific zenith angle Z corresponding to lati- 
tude <f>(I,J,t) the absorbed insolations A(2,6,t) and A(6,10,t) 
of Eqs . (4-28) and (4-29) have been determined for the given 

time t of 16 October 1973 for which the specific hour angles 
were applicable. The instantaneous parameter REF(t) of Eq . 
(4-35) was determined by the equation 
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REF(t) = F(2 ) - A26 - A610 - ( IA10 + ISIO) 



(6-2) 



where the expression in parentheses is identical to the in- 
stantaneous absorbed insolation of earth given by the numer- 
ator of ( 4-37 ) . 

The net terrestrial flux F 2 * was derived in Eq. (3-11) 
and was not subject to any systematic diurnal variation. 

Thus F 2 * may be taken as constant for a 24-hour period. 

Note that the first symbol at k=2 in Fig. 6 is not F(2), 
but QAVE, defined for each gridpoint by the analog of (4-1), 
namely 

QAVE = S(— )~ 2 Cos Z (6-3) 

r m 

Here Cos Z is the 24-hour average of Cos Z, and is given by 
Cos Z = [H Sin <j> Sin 6 + Cos <j> Cos 6 Sin H] / tt (6-4) 

where 

H = arc Cos [-Tan <f> Tan 6] (6-5) 

is the hour angle from local noon to sunset at the gridpoint. 
Thus the instantaneous values of incident insolation and of 
reflected insolation have been converted from F(2) and REF 
to their 24-hour averages denoted by 

QAVE = F(2) Cos Z/Cos Z (6-6) 

and 

QREF = REF Cos Z/Cos Z (6-7) 

respectively. 
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The reason for the conversion in (6-6) and (6-7) is to 
convert the net incoming insolation at k=2 to values that 
are comparable in time-scale to those of F 2 * . This device 
has the advantage of showing whether the resultant averaged 
heat-transf er quantities are compatible with the expected 
approximation to a heat budget for the given date. However, 
if it were desired to exhibit the expected diurnal variation 
of heat-transfers across the layers (2,6), (6,10) and at the 
ground, the quantity QAVE and all other solar-heating rates 
to be discussed would be multiplied by the factor 

(Cos Z/Cos Z) 

assuming that all the layer absorptivit ies and the solar re- 
flectivity are constant during the day for which the diurnal 
variation of heat transfers are of interest. 

With the solar quantities averaged for the day, one ob- 
tains 

BALT = QAVE - (QREF + F 2 *) (6-8) 

as a 24-hour average radiational balance in cal cm -2 min -1 
for the tropopause. 

C. CROSS-SECTIONAL RADIATIVE TRANSFERS IN LAYER (2,6) 

At this layer, in Fig. 6, the heat-transfers are assumed, 
as in the FNWC heating model, to be radiative only. For con- 
venient comparison with climatology, the heating (cooling) 
rate is given by 

BAL26 = Q26 - F26 (6-9) 
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where Q26 is the daily average absorption in layer (2,6) and 
is defined relative to A26 by an equation analogous to Eq. 
(6-6). The terrestrial cooling rate F26 is defined by Eq. 

( 3-12a) . 

Note that all values of BAL26 were negative in all cross- 
sectional diagrams, Figs. 7,..., 10. 

D. CROSS-SECTIONAL DATA FOR LAYER (6,10) 

1 . Radiative Transfers Only 

Neglecting temporarily the turbulent transfers at 
the earth-air interface, the present heating model computes 
the 24-hour average radiative cooling computed by the differ- 
ence equation 

BAL610 = Q610 - F610 (6-10) 

in a manner analogous to Eq. (6-9). The parameter BAL610 
was negative at all gridpoints of Figs. 7,..., 10. In fact 
the ratio R = F610/Q610 varied with latitude as follows 

0° 10° 20° 30° 40° 50° 60° 

R = 1.916 2.219 2.577 2.893 3.912 4.122 4.903 

Even when the instantaneous (near-noon) radiation balance 
was considered, the difference remained at least slightly 
negative in low latitudes corresponding to cooling. In high- 
er latitudes, instantaneous cooling rates were never reversed 
by consideration of the instantaneous computation of the ra- 
diation balance in layer (6,10). 
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2. Heat Balance at the Earth’s Surface 



In Section V.B, the earth-to-air turbulent transfer 
of latent heat, E, has been modeled after Kaitala (1974), 
through Eq. (5-1). The sensible heat transfer, Hp , given by 

H p = rE 

where r is the Bowen ratio determined empirically by Eq. 

(5-5) as a function of latitude. 

The sum of the turbulent heat transfers, Hp+E, has 
been subtracted from the average radiation at the earth's 
surface 

R = QABG - Fj o* (6-11) 

to form the average warming rate, BALB, at the surface 

BALB = QABG - F 10 * - (Hp+E) (6-12) 

where QABG is related to ABG of Eq. (4-37) by an analog to 
Eq. (6-6) and Fi 0 * is defined by Eq. (3-10) which is taken 
as constant for a 24-hour period as was the case with F 2 *. 

The parameter BALB has been shown in each cross-sec- 
tion as a function first of gridpoint on each meridian, and 
subsequently interpolated to 5° latitude increments as indi- 
cated in connection with Eq . (6-1). 

3 . Sensible Heat Modification of the Atmospheric Columns 

In Figs. 7,..., 10 the sensible heat transfer, Hp , 
which crosses the earth-air interface has been added to the 
air only in the layer (8,10), (cf. Kaitala, 1974). In this 
connection, the net radiative loss in the layer (6,10) has 
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been subdivided equally into two halves, BAL68 and BAL810 
(excluding Hp ) , equal to half the radiative loss BAL610. 

The sensible heat transfer Up , which has been allowed 
to diffuse only into layer (8,10), is shown as item (n) of 
layer (8,10) in Fig. 6. In the numerical computations, Hp 
is the smallest of all the heat-transfers considered in any 
atmospheric column presented in Figs. 7,..., 10. 

E. MERIDIONAL CROSS-SECTIONS OF THE VERTICAL HEAT BUDGET 

Note that the single-time set of heating computations 
for each of the four meridians, 125°W, 170°W, 145°E, and 
35°W, are shown in Figs. 7,..., 10 respectively. These fig- 
ures show the individual computations of the warming and 
cooling which are schematically formatted in Fig. 6, together 
with the layers to which they contribute in the heating mod- 
el . 

Note also that the four meridional cross-sections de- 
picted each have two sections: section (a) depicts results 
for a "tropical" section with gridpoints extending from 20°S 
(or the southernmost latitude to which the meridian in ques- 
tion extends) to 25°N, while section (b) of Figs. 7,..., 10 
depicts the gridpoint computations (averaged for 24-hours) 
over the latitudinal range 30°N to 65°N (or the northernmost 
latitude to which the meridian in question extends). Recall 
the latitudinal boundaries for the selected meridians were 
required to remain over the ice-free ocean on the FNWC grid. 

The latitudinal variation of both the short- and long- 
wave calculations, as observed in Figs. 7,..., 10, seem to 
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have uniform gradients across each meridian in spite of the 
differing climatic zones represented. Hence it was reason- 



able to 
October 



formulate a mean radiational balance model for 16 
1973, which is presented in Section VII. 
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LAT. latitude of the gridpoints in 







X i— I 
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Figure 6. Key to Meridional cross-sections for Figs. 7,..., 10. All radiative and 
turbulent transfer values are in cal cm -2 min -1 for the levels or layers considered. 
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Fig. 6 for key. Values computed from data for 16 October 1973. 
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October 1973. 



6206 0.6078 0.6005 0.5816 0.5533 0.5309 0.4996 
2816 0.3343 0.3430 0.3089 0.2861 0.2677 0.2202 
3099 0.2959 0.2992 0.3525 C.3963 0.4004 0.3943 
0292 -0.0224 -0.0417 -0.0793 -0.1241 -0.1373 -0.1150 
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Figure 9(a). 145°E Longitudinal cross-section, tropical section. Refer to Fig. 

for key. Values computed from data for 16 October 1973. 
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key. Values computed from data for 15 October 1973. 
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VII. MEAN BUDGETARY BALANCE MODEL 



By averaging key parameters, discussed in' Section VI, 
across the four possible meridional lines at each 5° latitude 
increment from 0° to 65° N, a mean heat budget model was con- 
structed which may be considered representative of Fall in 
the annual cycle. The advantage of such averaging is to fa- 
cilitate the comparisons of this study's results with other 
published heat budget results (cf. Budyko, 1956) even though 
such results are generally annual means. Also, by averaging 
over the four meridians studied here, meridional differences 
tend to be smoothed out, so that the zonally-averaged heat 
budgets of this study may be compared with others in the lit- 
erature . 

A. RADIATION MODEL OF THE OCEAN-ATMOSPHERE SYSTEM 

Fig. 11(a) depicts the schematic format of the mean ra- 
diation meridional cross-section (Figs. 11(b) and (c)). 
Turbulent transfer parameters have been omitted from this 
display and only Northern Hemisphere values are given. 

Note that the radiational quantities for the layers 
(6,8) and (8,10) of Figs. 7,..., 10 have been combined at 
identical latitudes in Fig. 11 as for the layer (6,10) and 
that all resultant values are meridionally averaged where 
possible. Since only the 170°W meridian contains results at 
65°N, the value at this latitude is not representative of the 
zonal-averaging process described above. 
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LAT latitude of the gridpoints in 5° increments 
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All radiational values have been averaged over a 24-hour day and are expressed in 

cal cm -2 min -1 for the levels or layers considered. The over-bar symbol refers to mean 

values averaged over the possible meridional lines (125°W, 170°W, 145°E, and 35°W) . 
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data for 16 October 1973. 







^ X H 



O 

P 

u 

0 

0 

ct: 



w 

0 

T) 

G 

P 

•H 

P 

G 

t-H 

I 

T3 

■H 

6 

Jh 

O 



G 

O 

♦H 

P 

0 
0 
CO 

1 

w 

CO 

o 

u 

o 



ci 

G 

O 



U 

0 

6 

c 

o 

•H 

P 

G 

*H 

G 

U 



G 

G 

0 



O 



0 

Jh 

G 

b£> 

•H 

Ip 



88 



Fig. 11(a) for key. All values listed are in cal cm -2 min -1 and are computed from 
data for 16 October 1973. 



To facilitate a graphical display of the results con- 
tained in Fig. 11 and for ease of comparison with published 
yearly climatological results (Malkus, 1962), the model dis- 
played in Fig. 11 has been used as the basis for testing the 
following component heat balances: (a) the earth-atmosphere 
system balance (using symbols consistent with Malkus (1962)), 
(b) the atmospheric radiative balance, and (c) the earth- 
surface balance. 

B. RADIATION BALANCES OF THE OCEAN-TROPOSPHERE SYSTEM 

The first graph, Fig. 12, consists of a radiation bal- 
ance composed of three parameters, R g , R , and R, which are 

defined (after Malkus, 1962) as follows: (1) R is mean ra- 

s 

diative energy gain (loss) rate of the ocean-troposphere 
system which was referred to as BALT at k=2 in Fig. 11(a); 

(2) R is the mean radiative cooling rate in the troposphere 

cl 

which is defined as the sum of BAL26 and BAL610 from Fig. 
11(a); and (3) R, the mean radiative warming (cooling) rate 
at the earth's surface which was identified as BAL10 in Fig. 
11(a). The relationship between the three radiative param- 
eters can be expressed as 

R = R + R (7-1) 

s a 

Table III lists quantitative comparisons of R , R , and 

S cl 

R, as found in this study for 16 October 1973, with values 
reported in Malkus (1962) for the climatological year, and 
values reported by Jenks (1974) for 25 April 1973. The lat- 
itudinal ly-averaged values are obtained by weighting each 
value by the cosine of the latitude 
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Figure 12. Radiational balance at the tropopause (R ); in 

o 

the atmospheric column between k=2 and k=10 (R ) ; and at the 

cl 

surface (R) for 16 October 1973. 
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(a) 


(b) 


(c) 


°LAT 


R / R / R 

S 1 S 2 S 3 


R / R / R 

a i a 2 / a 3 


/ ^2 / ^3 



0 


.0456/. 0627/. 031 


-.1231/-. 1235/-. 132 


. 1687/ 


/ . 154 


5 


. 0358/ . 0608/ . 024 


-.1361/-. 1254/-. 142 


.1720/. 2033/. 158 


10 


.0149/. 0589/. 032 


-.1463/-. 1292/-. 167 


. 1596/ 


/ . 194 


15 


-.0060/. 0513/. 041 


-.1522/-. 1387/-. 167 


.1463/. 2090/. 208 


20 


-.0353/ . 0437/. 032 


-.1713/-. 1368/-. 179 


. 1361/ 


/ . 206 


25 


-.0620/. 0361/. 062 


-.1896/-. 1235/-. 174 


. 1276/ . 1900/ . 236 


30 


-.0871/. 0266/. 068 


-.1853/-. 1178/-. 160 


.0982/ 


/. 225 


35 


-.0944/. 0171/. 055 


-.1650/-. 1102/-. 153 


.0705/. 1482/. 201 


40 


-.1066/ 0.0 / . 044 


-.1770/-. 1159/-. 140 


.0703/ 


/ . 175 


45 


- . 1240/- .04 18/ -.003 


-.1738/-. 1273/-. 135 


.0499/. 1026/. 132 


50 


-.1423/- .0570/. 009 


-. 1759/-. 1330/-. 128 


. 0336/ 


/ . 124 


55 


- . 1325/- .0684/ . 022 


-.1472/-. 1330/-. 132 


.0146/. 0646/. 155 


60 


-.1517/- .0760/. 042 


-.1528/-. 1311/-. 141 


.0013/ 


/ . 184 


65N 


-.1112/ / . 106 


-.1233/ /- . 114 


. 0120/ 


/ . 219 


Wt. 

Avg. 


- . 0552/ . 0186/ . 039 


-.1589/-. 1264/-. 150 


.1037/. 1625/. 185 



Table III. Comparisons of Radiational Balances (a) at the 
tropopause R , (b) in the vertical column R between k=2 and 

S ci 

k=10, and (c) at the earth's surface R of the meridionally 
averaged values as found in this study denoted by subscript 
"1". The corresponding annual values reported by Malkus 
(1962) are denoted by subscript "2", and those for 25 April 
1973 found by Jenks (1974) denoted by subscript "3". All 
values are in ly min -1 . 
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Wt. Avg. 



(7-2) 



n 

Q . Cos d> . 

i=l c i _i 

n 

. Cos d> . 

1=1 y l 

It can be noted from Table III that although the R, val- 

s 

ues found by this study are consistently lower than those re- 
ported by Malkus for the year, the seasonally higher values 
reported by Jenks tend to offset this difference. The same 
holds true for the balance at the surface (R). However, 
both seasonal studies show greater radiative cooling in the 
troposphere (R ) than the yearly climatological values which 

ci 

evidently was obtained after Budyko (1956) as a residual be- 
tween R g and R, as functions of latitude. It is noteworthy 

that Jenks' (1974) R results (column (c) of Tabic III) are 

s 

greater than those obtained here by .094 ly min -1 (in the 
mean), and of this surplus the additional heating rate 
AR = 0.081 ly min -1 is found at the surface. 

C. TROPOSPHERIC HEATING 

Fig. 13 graphically depicts the mean tropospheric heat 
budget as found in this study. The parameters presented in- 
clude: R which was discussed in Sec. VI I. B, (E+H r ) which is 
the turbulent heat transfer rate due to latent and sensible 
heat transport realized in the atmosphere, Q va which is the 
net heat flux divergence rate in the troposphere required to 
bring about a heat balance (which would not be expected to 
be achieved for the time scale here), and S the storage 

<X 

term. These terms are related by 
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heat loss in the Troposphere (Q va +S ) is depicted as the sum 

of the tropospheric radiative cooling (R ) and the turbulent 

a 

transfer heating (E+Hp) in the troposphere. 
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LAT 


(a) 

R / R 

a i a 2 


(b) 

(E+H r ) i /(E+H r ) 2 


(c) 

(Q +S ) / (Q +S ) 
v va a i ' v va a z 


0 


-.1231/-. 1235 


.0294/ 


-.0937/ .0874 


5 


-. 1361/-. 1254 


.0332/ .1539 


0.1029/ .0285 


10 


-.1463/-. 1297 


. 0425/ 


-.1038/ .0437 


15 


-.1522/-. 1387 


. 0636/ . 1786 


-.0916/ .0095 


20 


-.1713/-. 1368 


.0544/ 


-.1169/-. 0057 


25 


-. 1896/-. 1235 


.0744/. 1748 


-.1095/- . 0095 


30 


-.1853/-. 1178 


. 1164/ 


-.0749/ .0190 


35 


-. 1650/-. 1102 


. 1701/. 1615 


-.0051/ .0190 


40 


-.1770/-. 1159 


.1544/ 


-.0226/ .0038 


45 


-.1738/-. 1273 


.1518/. 1254 


-.0220/-. 0095 


50 


-.1759/-. 1330 


.1486/ 


-.0273/-. 0152 


55 


-. 1472/-. 1330 


.0880/. 0988 


-.0592/-. 0228 


60 


-.1528/-. 1311 


. 0924/ 


-.0604/-. 0266 


65N 


-.1233/ 


. 0188/ 


- . 1045/ 



wt 

Avg. 1589/-. 1264 . 0855/ . 1533 -.0735/ .0201 

Table IV. Comparisons of trapospheric radiative cooling rate 

(R ), turbulent transport warming rate (E+H ) , and net trop- 
a r 

ospheric heat loss rate (Q +S ) as found by this study for 

va a 

16 October 1973 denoted by subscript "1" and as reported by 
Malkus (1962) for yearly climatology denoted by subscript "2". 
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Q + S = R + (E+H„ ) 
va a a v r 



(7-3) 



The storage rate S a has not been separated specifically from 
Q va , but over the annual cycle it is assumed zero. Assuming, 
however that Q ya = 0> if the weighted average values from 
Table IV are set into Eq. (7-3), to determine the mean stor- 
age cooling rate, S , the latter turns out to be 

cl 



S a = (|^) = - . 54° C (day)" 1 

for the time (16 October 1973) under study. The net lateral 
flux-divergence term Q va has not been assessed here, but pre- 
sumably its mean effect would show up as cold advection 
across the equator, and would result in a reduced value of 




In Fig. 13, the dashed portion of the (E+IIp) and (Q va + ^ a ) 
curves are displayed to show the effects of varying the K* 
value in the calculation of latent heat transport (Eq. 5-1), 
as discussed in Section V. In the dashed portion, the K* 
value of 2 x lo 5 cm 2 sec -1 was tested for latitudes 0°, 5°, 
and 10°N, with a transitional value of 1.5 x io 5 cm 2 sec -1 
used at 15°N. At all points on the solid (E+Hp) line K* was 
1 x io 5 cm 2 sec -1 . This smaller K* value seemed to provide 
more realistic values of (E+Hp) in mid-latitudes but to some- 
what unrealistically small values in low latitudes. 

A possible explanation for the small (E+Ilp) values at 
low latitudes is the difficulty in correctly determining the 
temperature and humidity gradients between the surface and 
the top of the constant flux layer in the data-sparse 
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equatorial areas. In these tropical latitudes some time- 
and space-smoothing of the few available radiosondes could 
have been introduced into the original FNWC map analyses, 
which would tend to smooth out vertical gradients as against 
those to be expected with observed radiosondes on 16 October 
1973. Hence it was considered justifiable to revert to 
Kaitala's (1974) suggested K* = 2 x 10 s cm 2 sec -1 in low 
latitudes, 0°-10°N, where the resultant dashed (E+H^,) curve 
of Fig. 13 qualitatively appears more reasonable. For ex- 
ample, the peak in (E+H^,) at 15°N appears to be due to rela- 
tively strong evaporative and sensible heat transfers 
associated with the Trade Wind regime, while the higher 
latitude peak (35°-50°N) in Fig. 13 is due to the trend to- 
wards winter maximum turbulent-heat transfers associated 
with cool continental air flowing over warm ocean currents. 

No attempt to select a final K* value is suggested from 
these results. If the same turbulent-transfer model can be 
tested over a large enough set of seasonal samples, it should 
be possible to assess the appropriate K*-value which yields 
the best approximation to the atmospheric heat balance. 

D. SURFACE HEATING RATES 

Fig. 14 ""depicts the zonally-averaged surface heating bud- 
get as a function of latitude resulting from this study. 

The parameters include the net surface radiative warming (R), 
the surface heat loss due to turbulent transfer, -(E+H^), 
and the net warming or cooling at the surface, (Q +S). 

These parameters are related by the following expression 
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20 




-20 i i i i i i i i i i i i i 

0 5 10 15 20 25 30 35 40 45 50 55 60 65 

Latitude ( N) 

Figure 14. Surface heat budget disposition. The net warming 
(cooling) rate at the earth's surface (Q^.^+S) is given by the 
sum of the net radiative surface warming rate R and the loss 
rate due to turbulent heat transfer, -(E+H). 
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LAT 


(a) 

R, / R 2 


(b) 

-(E+Hp) i /-(E+H r ) 2 


(c) 

(Q +S) / (Q +S) 

VO l VO 7 2 


0 


. 1687/ 


-.0294/ 


.1393/ 


5 


.1720/. 2033 


-.0332/-. 1539 


.1388/ .0494 


10 


. 1596/ 


- . 0425/ 


.1171/ 


15 


. 1463/. 2040 


-.0636/-. 1786 


.0827/ .0304 


20 


.1361/ 


-.0544/ 


.0817/ 


25 


.1276/. 1900 


-.0774/-. 1748 


.0502/ .0152 


30 


. 0982/ 


-.1104/ 


-.0122/ 


35 


.0705/. 1482 


-.1701/-. 1615 


-.0996/-. 0133 


40 


. 0703/ 


-.1544/ 


-.0841/ 


45 


.0499/. 1026 


-.1518/-. 1254 


-.1019/-. 0228 


50 


. 0335/ 


-.1486/ 


-.1151/ 


55 


. 0146/ . 0646 


-.0880/-. 0988 


-.0734/-. 0342 


60 


. 0013/ 


-.0924/ 


-.0911/ 


65N 


. 0120/ 


-.0188/ 


- . 0068/ 



wt . 

Avg. .1037/. 1625 0855/- . 1533 .0182/ .0092 

Table V. Comparisons of net radiative warming rate at the 
surface (R), the heat loss rate at the surface due to turbu- 
lent transport -(E+Hp), and the net warming or cooling at the 
earth's surface (Q +S) for 16 October 1973 (denoted by sub- 
script "1"). For comparison, annual mean values reported by 
Malkus (1962) are listed at selected latitudes (denoted by 
subscript "2"). All values are in ly min -1 . 
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Q vo + S = R - (E+H r ) (7-4) 

where S is the oceanic heat-storage term. Again it was not 
possible here to separate the contributions, Q vq and S. The 
square points and dashed-in portion of the -(E+II^) and 
curves depict the effect of varying the K* value in Eq. (5-1) 
as discussed in Sec. VII. C. The overall distribution in the 
net oceanic heating rate (Q vq +S) is the same for both choices 
of K* . However this rate is somewhat reduced for the case 
of the larger K* values suggested as applicable for the 
tropics . 

Table V shows quantitatively the difference in the param- 
eters depicted in Fig. 14 as found by this study for 16 
October 1973 and the climatological annual values reported 
by Malkus (1962). Note the smaller mean seasonal radiation- 
al warming rate (R) and the effect of the smaller mean heat 
loss by turbulent transfer in giving a slightly larger mean 
warming of the surface. 
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VIII. CONCLUSIONS 



A heat budget model for inclusion in the FNWC numerical 
weather prediction scheme was evaluated using data along 
four meridional lines over the Pacific and Atlantic Oceans. 
These data soundings were used to compute a space-time aver- 
age of the radiational and turbulent heating parameters for 
one day over an ocean environment. This space-time average 
agreed reasonably well with values derived from other stud- 
ies, including satellite measurements. There were some ob- 
vious differences which could be ascribed to the seasonality 
of the data used in this study comprising only a one day 
space-time average while the comparison data was averaged 
over years to arrive at climatological results. 

Other differences could be ascribed to the two-cloud 
layer parameterization for cloud amounts, the treatment of 
cloud reflectivity and absorptivity for solar wave lengths, 
and the parameterization of sensible and latent heat trans- 
fers. At this stage in the model testing, it would seem 
that the radiational estimates are still somewhat question- 
able due to the adoption of the Smagorinsky-type cloud 
amounts. However the regular variation of the radiative 
terms from gridpoint to gridpoint tends to support the ap- 
parent realism of the radiative model. The turbulent trans- 
fer heating rates of the model seem to be subject to a 
greater degree of possible random error due to uncertainty 
of the vertical gradients involved near the surface. 
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To complete the comparisons of the radiational aspects 
of this model with cited climatological data (see Sec. VII) 
and with satellite data, further testing should be done at 
different seasons of the year and possibly over other oceanic 
areas. Thus with a full year of representative data, refine- 
ments of the cloud-cover parameterization and in the cloud 
absorption and reflection values used in the model may be 
assessed . 

Finally, the greatest gridpoint to gridpoint variability 
occurred in connection with the computation of evaporative 
and sensible heat parameters. It is perhaps here that the 
moisture and heating source terms used in the FNWC predic- 
tion model might best be improved. 
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